arXiv:1508.06920v2 [math.NT] 28 Aug 2015 


DESINGULARIZATION OF COMPLEX MULTIPLE ZETA-FUNCTIONS 


HIDEKAZU FURUSHO, YASUSHI KOMORI, KOHJI MATSUMOTO, AND HIROFUMI TSUMURA 


Abstract. We introduce the method of desingularization of multi-variable multiple zeta- 
functions (of the generalized Euler-Zagier type), under the motivation of hnding suitable 
rigorous meaning of the values of multiple zeta-functions at non-positive integer points. We 
reveal that multiple zeta-functions (which are known to be meromorphic in the whole space 
with infinitely many singular hyperplanes) turn to be entire on the whole space after taking 
the desingularization. The desingularized function is given by a suitable finite ‘linear’ com¬ 
bination of multiple zeta-functions with some arguments shifted. It is shown that specific 
combinations of Bernoulli numbers attain the special values at their non-positive integers 
of the desingularized ones. We also discuss twisted multiple zeta-functions, which can be 
continued to entire functions, and their special values at non-positive integer points can be 
explicitly calculated. 


0. Introduction 


We begin with the multiple zeta-function of the generalized Euler-Zagier type 

defined by 

OO oo r 

Cr((si);(7i)) = Cr(si,---,Sr;7i,--->7r) := Y[ -^ ( 0 - 1 ) 

mi=l 771^ = 1 ^ = 1 

for complex variables si,... ,Sr, where 71 ,... , 7 ^ are complex parameters whose real parts 
are all positive. Series ( |0.1[ ) converges absolutely in the region 

T>r = {(si, . . . , Sr) G I 3^(Sr-fc-|-l + ■ ■ ■ + Sr) > k {1 ^ k ^ r)}. (0.2) 


The first work which established the meromorphic continuation of (0.1) is Essouabri’s thesis 
[5]. The third-named author [HI Theorem 1] showed that ( |0.1[ ) can be continued meromor- 
phically to the whole complex space with infinitely many (possible) singular hyperplanes. 

A special case of (0.1) is the multiple zeta-function of Euler-Zagier type defined by 


OO r 

Cr((Si)) =Cr(si,S2,---,Sr) = ^ H -^ "ij)”*" > (0-3) 

mi,...,mr=l j=l 

which is absolutely convergent in 

Note that Cr(('Sj)) = Cr(('Sj); (1)). Its special value Cr{ni,... ,nr) when ni,... ,nr are pos¬ 
itive integers makes sense when > 1. It is called the multiple zeta value (abbreviated as 
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MZV), history of whose study goes back to the work of Euler [ 6 ] published in 1776 For 
a couple of these decades, it has been intensively studied in various fields including number 
theory, algebraic geometry, low dimensional topology and mathematical physics. 

On the other hand, in the late 1990s, several authors investigated its analytic properties, 
though their results have not been published (for the details, see the survey article [T 8 ]h 
In the early 2000s, Zhao [25] and Akiyama, Egami and Tanigawa [T| independently showed 


that (0.3) can be meromorphically continued to C^. Furthermore, the exact locations of 


singularities of (0.3) were explicitly determined in [Tj: Cr{{sj)) for r > 2 has infinitely many 
singular hyperplanes 

Sf — 1, Sj.—! + Sj- — 2,1,0, 2, 4, 6 ,..., 

Sr-k+i + Sr-k+2 H- 'c Sr = k - Tl {3 ^ k ^ r, n ^ No). (0.4) 

It is natural to ask how is the behavior of Cr(~R'i) ■ ■ ■, —nr) when ni,..., are positive (or 
non-negative) integers. However, unfortunately, almost all non-positive integer points he on 
the above singular hyperplanes, so they are points of indeterminacy. For example, according 

to [Il| 2 |, 

lim lim C 2 (ei,e 2 ) = 

S \—>^0 £2—^0 O 

5 

lim lim C 2 (£ 1 , 62 ) = —, 

£2 —^0£i—>^0 Iz 

3 

lirnC 2 (£,e) = -• 

e-S’O 8 

There are some other explicit formulas for the values at those non-positive integer points as 
limit values when the way of approaching to those points are fixed ([Il[2l|22l|23l[l5l|2l]). 

However points of indeterminacy cannot be easily investigated, so we can raise the following 
fundamental problem. 

Problem 0.1. Is there any ‘rigorous’ way to give a meaning of Cri—ni ,..., —Ur), without 
ambiguity of indeterminacy, for m,..., G 

Several approaches to this problem have been done so far. Guo and Zhang m, Manchon 
and Paycha m and also Guo, Paycha and Zhang |10j discussed a kind of renormalization 
method. In the present paper we will develop yet another approach, called the desingulariza- 
tion, in Section]^ The Riemann zeta-function (^{s) is a meromorphic function on the complex 
plane C with a simple and unique pole at s = 1. Hence (s — l)C('S) is an entire function. This 
simple fact may be regarded as a technique to resolve a singularity of C('S) and yield an entire 
function. Our desingularization method is motivated by this simple observation. For r ^ 2, 
multiple zeta-functions have infinitely many singular loci. We will show that a suitable finite 
sum of multiple zeta-functions will cause cancellations of all of those singularities to produce 
an entire function whose special values at non-positive integers are described explicitly in 
terms of Bernoulli numbers (see Figure and (4.3) for the case r = 2). 


^ You can find several literatures which cite the paper saying as if it were published in 1775. But according 
to Euler archive http://eulerarchive.maa.org/, it was written in 1771, presented in 1775 and published in 
1776. 
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Singularities of 

C2('Si,S2) - 

C2(Si — 1, S2 + 1) - 

^2(51 — 2 , S2 + 2 ) - 


Figure 1. Singularities of C 2 ’s 


Another possible approach to the above Problem 0.1 is to consider the twisted multiple 


series. Let ^ 1 ,..., G C be roots of unity. For 71 ,..., 7 ^ G C with 5 f 7 j > 0 (1 ^ j ^ r), 
define the multiple zeta-function of the generalized Euler-Zagier-Lerch type by 

00 CO r 

Cr((sj);(0);(7i)) := + + (0.5) 

mi=l mr=l j=l 


which is absolutely convergent in the region 2?^ defined by (0.2). We note that the multiple 


zeta-function of the generalized Euler-Zagier type (0.1) is its special case, that is, 


Cr{{sj); (7i)) = Cr{{sj)-, (1); (7i))- 


Because of the existence of the twisting factor ^ 1 ,..., we can see (in Theorem 2.1 below) 


that, if no is equal to 1, series (0.5) can be continued to an entire function, hence its values 


at non-positive integer points have a rigorous meaning. Moreover we will show that those 
values can be written explicitly in terms of twisted multiple Bernoulli numbers. 

In Section we will introduce multiple twisted Bernoulli numbers, which are connected 


with multiple zeta-functions of the generalized Euler-Zagier-Lerch type (0.5). After discussing 


the aforementioned properties of (0.5) in Section]^ we will develop our method of desingu- 
larization in Section Multiple zeta-functions ( |0.1[ ) are meromorphically continued to the 
whole space with their singularities lying on infinitely many hyperplanes. Our desingular- 


ization is a method to reduce them into entire functions (Theorem 3.4). We will further 


show that the desingularized functions are given by a suitable finite ‘linear’ combination of 


multiple zeta-functions (0.1) with some arguments shifted (Theorem 3.8) This is the most 
important result in the present paper, in which we see a miraculous cancellation of all of their 
infinitely many singular hyperplanes occurring there by taking a suitable finite combination 
of these functions. We will also prove that certain combinations of Bernoulli numbers attain 


the special values at their non-positive integers of the desingularized functions (Theorem 3.7). 
Several explicit examples of desingularization will be given in Section]^ 

It is to be noted that these observations on our desingularization method lead to the 
construction of p-adic multiple L-functions which will be discussed in a separate paper [ 8 ] . 
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1. Twisted multiple Bernoulli numbers 


In this section, we first review the definition of classical Bernoulli numbers and Koblitz’ 
twisted Bernoulli numbers. Then we will introduce twisted multiple Bernoulli numbers, 
their multiple analogue, and investigate their expression as combinations of twisted Bernoulli 
numbers. 

Let N, No, Z, Q, M and C be the set of natural numbers, non-negative integers, rational 
integers, rational numbers, real numbers and complex numbers, respectively. For s G C, 
denote by JRs and Ss the real and the imaginary parts of s, respectively. 

It is well-known that ((s) is a meromorphic function on C with a simple pole at s = 1, 
and satisfies 



(k€N>i) 
(k = l), 


( 1 . 1 ) 


where {B^} are the Bernoulli numbers defined by 



(see [211 Theorem 4.2]). 

Definition 1.1 f |131 p.456]). For any root of unity we define the twisted Bernoulli 
numbers {iBn(^)} by 



n=—l 


where we formally let (—1)! = 1. 


Remark 1.2. Koblitz m generally defined the twisted Bernoulli numbers associated with 


primitive Dirichlet characters. The above {®n(0} correspond to the trivial character. 
In the case .^ = 1, we have 


®_i(l) = -l, 5S„(l) = -i^ (neNo). 

n -|- 1 


(1.3) 


In the case ^ / 1, we have iB_i(^) = 0 and 53^(0 = ^ where {Hn{X)}n^o 

are what is called the Frobenius-Euler numbers associated with A defined by 



n=0 


(see Frobenius El)- We obtain from (|1.2[) that ^n{C) ^ Q(0- For example. 



(1.4) 


^ or better to be called Seki-Bernoulli numbers, because Takakazu (Kowa) Seki published the work on these 
numbers, independently, before Jakob Bernoulli. 
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Let Hk be the group of A:th roots of unity. Using the relation 


1 


k 




1 


X -1 - 1 ^ 1 - 


{k G N>i) 


(1.5) 


for an indeterminate X, we obtain the following. 

Proposition 1.3. Let c G N>i. For n G No, 

Bn+l 


(1 - ^ = E *»(«)■ 


( 1 . 6 ) 


5^=1 

5#i 


Remark 1.4. Let ^ be a root of unity. As an analogue of (1.1), it holds that 


cP{-k-0 = ^k{0 (feGNo), (1.7) 

where (j)(s; ^) is the zeta-function of Lerch type dehned by the meromorphic continuation 
of the series 


m^l 


(cf. [m Chapter 2, Section 1]). 


We see that (1.6) can also be given from the relation 

(c^-^ -1) c(s) = Y1 


(1.9) 


5^=1 

5/1 


Now we define certain multiple analogues of twisted Bernoulli numbers. 

Definition 1.5. Let r G N, 71 ,..., 7 ^ G C and let ^ 1 ,..., G C \ {1} be roots of unity. Set 

( 1 . 10 ) 


1 


^r-((ti);(^i);(7i)) := = n — 

i=i k=j j=i 1 “ exp ^ 7 j Ylk=j 

and dehne twisted multiple Bernoulli numbers {*B(ni,.. ., (^j); ( 7 j))} by 

j-Fll >nr 

(fl); (7,)) = i8(ni,..., n,; (f,); (7,))-!- ■ ■ ■ E... (1.11) 

Til'. ^r' 

ni={j nr=0 

Remark 1.6. It is possible to generalize the above definition to the case when = 1. In this 
case, the sum with respect to on the right-hand side of ( 1 . 11 ) is from —1 to oo, hence 
gives a more natural extension of (|1.2[). 


In the case r = I, we have IBn(^i) = lB(n;^i; I). Note that since / 1 (1 ^ j ^ r), we 
see that Sjr{{tj); {^j); {'Jj)) is holomorphic around the origin with respect to the parameters 
ti,... ,tr, hence the singular part does not appear on the right-hand side of ( 1 . 11 ). 

We immediately obtain the following from (1.2), (1.10) and (1.11). 


'’We are not sure which is better, “twisted multiple”, or “multiple twisted”. But we will skip this problem 
because it looks that these two adjectives are “commutative” here. 
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Proposition 1.7. Let 71,...,7^ G C and G C \ {1} be roots of unity. Then 

®(ni,..., n,.; (^j); (7j)) can he expressed as a polynomial in { 5 Sn(^j) | 1 ^ J ^ f, n ^ 0} 
and {71,...,7r} with Q-coefficients, that is, a rational function in and {■jj} with Q- 
coeffcients. 


Example 1.8. We consider the case r = 2 . Substituting (1.2) into (1.10) in the case r = 2, 
we have 

-^2(ti,t2;6,6;7i)72) = 


1-6 exp (71 (ti + t2)) 1 - 6 exp (72t2) 


E ®m(6) 


iTiti + hY 


\rn=0 


ml 




/ 


Vn=0 


thi 


7£f 

n! 




fTl 

„n ''2 


m=0 n=0 

Putting I = n + j, we have 




k,j-^0 

k+j=m 


nl 


00 00 l 


^2(ti,t2;6,6;7i,72) = EEEC) 


k=Q 1 = 0 j =0 


kl ll ’ 


which gives 


6)71, 72 ) = E (0®fc+.(6)5Sz-i(6)7^'7E ^ No). (1.12) 

n \d / 


j =0 

For example, we can obtain from ( 1 . 4 ) that 

®(0,0;6,6;7d72) = 

®(o, 1 ; 6 ) 6 ; 71)72) = 


1 


(i- 6 )(i- 6 ) 

671 + 672 - 66(71 + 72) 


, *B( 1 , 0 ; 6 , 6 ; 71 ) 72 ) = 


671 


( 1 - 6 ) 2 ( 1 - 6 )’ 


®(i) 1 ; 6 ) 6 ; 71)72) = 


( 1 - 6 ) 2 ( 1 - 6)2 

?i 7 i( 7 i - 6(71+72)) +671(71 - 6(71 - 72)) 


( 1 - 6 ) 61 - 6)2 

The following series will be treated in our desingularization method in Section 
Definition 1.9. For c G M and 71,... ,7^ G C with iR-jj > 0 (1 ^ j ^ r), define 

^r((6); (7i);c) = n ' 


1=1 V exp 71 EL, 6-1 exp c-fj YIk=i 4-1 


/ 


r 00 


1=1 ^^=1 


1, EL, 4 ) 


m—1'' 


m\ 


(1.13) 


In particular when c G N>i, by use of ( | 1 . 5 [ ), we have 

r 

L(( 4 );( 7 i);c) = n E 


1=1 


1 -6exp 71 EL, 4 
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erT^l 

Remark 1.10. We note that S^riitj)] {'Jj)] c) is holomorphic around the origin with respect 
to the parameters (tj), and tends to 0 as c —)• 1. We also note that the Bernoulli numbers 
appear in the Maclaurin expansion of the limit 

These are important points in our arguments on desingularization methods developed in 
Section [3l 


Example 1.11. Similarly to Example 1.8, we obtain from (1.13) with any c G M that 


OO I I 


eIe(‘)(i 


k,l=0 i=o 


— c 


k+j+l 


1 - d' 


k+l 


l-j+i\ ^k+j+1 Bj-j+i ^k+j^i-j I 


k + j + ll - j + l 


7 i ^72 


k\l\ 


2 (1.15) 


Therefore it follows from (1.11) and (1.14) that 


^ ^ QS(/c,Z;6,6;7i,72) 


I 


j=0 

for c G N>i. 


E (') (l-(i-(Li. M„) 


(1.16) 


Remark 1.12. Kaneko m defined the poly-Bernoulli numbers {l?l^^}neNo ^ by use of 
the poly logarithm of order k. Explicit relations between twisted multiple Bernoulli numbers 
and poly-Bernoulli numbers are not clearly known. It is noted that, for example, 


= 


i 

E 

j=0 


j) J + 1 


which resembles ( 1 . 12 ) and (1.16). 


{I G No), 


2. Multiple zeta-functions 

Corresponding to the twisted multiple Bernoulli numbers {lB((nj); (^j); ( 7 j))} is the mul¬ 
tiple zeta-function of the generalized Euler-Zagier-Lerch type (0.5) defined in Introduction, 
which is a multiple analogue of This function can be continued analytically to the 

whole space and interpolates lB((nj); (.^j); ( 7 j)) at non-positive integers (Theorem 2.1). 

Assume 7 ^ 1 (1 ^ j ^ r). Using the well-known relation 

1 


u = 


r(s) 


e-ut^s-i 


dt, 


we obtain 

Criisj )] (^j); (7i)) 
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oo oo r / r 

E-E(n«r) n 



Hi 


r(sfc)y 7 ( 0 , 
ij^^v{-ij{TJk=jtk)) 


n exp(-4(X] n 


Sfe -1 


dti 




k=l 


J Ao,oor y 1 - Cj exp(-7i(ELi ^fc)) M 


n 


1 


n 


Cjexp(-7j(EL7^fc)) 


L^i (e 2 ™fc _ i)r(sfc) y 7e- 1 - exp(-7i(ELj ^fc)) 

e 




=(-ir n 


U (e2™fc _ i)r(sfc)^ 


^l»-((ti), ^), (7i)) n 


( 2 . 1 ) 


fc=i 


where C is the Hankel contour, that is, the path consisting of the positive real axis (top side), 
a circle around the origin of radius e (sufficiently small), and the positive real axis (bottom 
side). Note that the third equality holds because we can let e —)• 0 on the fourth member 
of ( |2.1[ ). In fact, the integrand of the fourth member is holomorphic around the origin with 
respect to the parameters (tj) because of 7 ^ 1 (1 ^ j ^ r). Here we can easily show that 
the integral on the last member of (2.1) is absolutely convergent in a usual manner with 
respect to the Hankel contour. Hence we obtain the following. 

Theorem 2.1. Let ... ,^r ^ C be roots of unity and 71 ,..., 7 ^ G C with iR'yj > 0 (1 ^ j ^ 
r). Assume that 

ij + 1 for all j (1 ^ j ^ r). ( 2 . 2 ) 

Then, with the above notation, Criisj)', (ij)) oan be analytically continued to as an 
entire function in (sj). For ni,..., G No, 


Cr((-nj);(Ci);(7j)) = (- 1 ) 


r+niH-hn. 


®(K);(^-^);(7i)). 


(2.3) 


Proof. Since the contour integral on the right-hand side of (2.1) is holomorphic for all (s^) G 
C^, we see that Cr-((sj); (ij)) can be meromorphically continued to and its possible 

singularities are located on hyperplanes G N (1 ^ A: ^ r) outside of the region 

of convergence because — l)r(sfc) does not vanish at Sk G Furthermore, for 

Sk = Ik ^ N, the integrand of the contour integral with respect to t^ on the last member of 
(2.1) is holomorphic around tk = 0. Therefore, for 1^ G N, we see that 


lirn / Siriitj), iCj ^), i'ljK ^dtk = I Sjri{tj), (^ ■ ^), (7i))Afc 'dtk = 0 , 

Sk^lk Je JCe 

because of the residue theorem, where = {ee*® | 0 ^ ^ 27r} for any sufficiently small 

e. Consequently this implies that Criisj)', ilj)) has no singularity on Sk = Ik, namely 
Cr(('Sj); (^j); ( 7 j)) is entire. Finally, substituting (1.11) into (2.1), setting (sj) = {—Uj) and 
using 


ih —1 . 


, 1 (-l)’"n! 

hm , „ . - , = - 

s^-n (e2’''*« — l)F(s) 27ri 


(n G No), 


we obtain (2.3). Thus we complete the proof of Theorem 2.1 


□ 


Such a type of explicit formulas for non-positive integer values of twisted multiple zeta- 
functions in several variables was already obtained by de Crisenoy [2] in a much more general 
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context (with real coefficients) by a quite different method. Some partial cases of Theorem 


2.1 are also recovered by the results in Matsumoto-Tanigawa m and Matsumoto-Tsumura 

m- 

In [T71 Theorem 1 ], it is shown that the multiple zeta-function Cr{{sj)', (Cj)] ilj)) of the 


generalized Euler-Zagier-Lerch type (0.5) with all = 1 is meromorphically continued to 


the whole space with possible singularities. A more general type of multiple zeta-function 


is treated in |15) . where equation (2.3) without the assertion of being an entire function is 


shown in the case of / 1 for all j and the meromorphic continuation of Cr((sj); {^j)] {ij)) is 
also given. We stress that in a separate paper [ 8 ], we construct a p-adic multiple L-function 
which can be regarded as a p-adic analogue of Cr{{sj)', (Cj)] {ij))- 


Remark 2.2. Without assumption (2.2), it should be noted that (2.1) does not hold generally, 


more strictly the third equality on the right-hand side does not hold because the Hankel 
contours necessarily cross the singularities of the integrand. 

In our recent paper j9], we will show the following necessary and sufficient condition that 
Cr{{sj)', (ij)) is entire, and will determine the exact locations of singularities when it is 

not entire: 

Theorem 2.3. Let ..., ^ C be roots of unity and 71 ,..., 7 ^ G C with > 0 (1 ^ j ^ 


r). Then Cr((sj); (fj); (jj)) can be entire if and only if the condition (2.2) holds. When it is 
not entire, one of the following cases occurs: 

(i) The function Criisj); {fj)', {'yj)) has infinitely many simple singular hyperplanes when 

= 1 for some j (1 ^ j ^ r — 1 ). 

(ii) The function Criisj)', (Cj)] (ij)) has a unique simple singular hyperplane = 1 when 

1 for all j (1 ^ j ^ r — 1 ) and = 1 . 


3. DESINGULARIZATION OF MULTIPLE ZETA-FUNCTIONS 

In this section we introduce and develop our method of desingularization. In our previous 
section we saw that the multiple zeta-function Criisj)', i'Jj)) of the generalized Euler-Zagier 


type ( 0 . 1 ) is meromorphically continued to the whole space with ‘true’ singularities whilst 


the multiple zeta function Criisj)', (^j); ( 7 j)) of the generalized Euler-Zagier-Lerch type (0.5) 


under the non-unity assumption ( 2 . 2 ) is analytically continued to as an entire function. 
Our desingularization is a technique to resolve all singularities of Cr(('Sj); ( 7 j)) to produce an 
entire function ( 7 j)). Consider the following expression; 

1 




(c _ 1)- 4- 

^1--L 


E 


Cr((si);(?i);(7j))- 


(3.1) 


This is surely nonsense, because c G N>i on the right-hand side. However, because of the 
holomorphy of Cr(('Sj); (?j); ( 7 i))j we observe that the left-hand side is also (at least formally) 
holomorphic. Our fundamental idea is symbolized in this primitive expression (3.1). Our 
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idea is motivated from a very simple observation 

(1 - s)C(s) = lini (c^"" - 1) C(s). 

c^l C — 1 

Here on the left-hand side we find an entire function (1 — s)((s), which is merely a product 
of (1 — s) and the meromorphic function ^(s) with a simple pole at s = 1. While on the right 
hand-side, when c G N>i, we may associate a decomposition 


into a sum of entire functions !)• 

Our desingularization method, a rigorous mathematical formulation to give a meaning of 
(3.1) will be settled in Definition 3.1 An application of desingularization to the Riemann zeta 

that our ( 7 j)) is 


3.4 


function (^(s) is given in Example 3.3 We will see in Theorem 
entire on the whole space C'’. We stress that ( 7 j)) is worthy of an important object 

from the viewpoint of the analytic theory of multiple zeta-functions. In fact, its values at not 
only all positive or all non-positive integer points but also arbitrary integer points are fully 
determined (see Example |4.9[ ). 

Theorem |3.7| will prove that suitable combinations of Bernoulli numbers attain the special 
values at non-positive integers of (jj))- Theorem 3.8 will reveal that our desingu- 

larized multiple zeta-function Cr^®((sj); (jj)) is actually given by a finite ‘linear’ combination 
of the multiple zeta-function Cr((sj + mj); (jj)) with some arguments appropriately shifted 


by rrij G Z (1 ^ j ^ r). Example |4.2| and Remark 4.3 are our specific observations for double 
variable case. 


Definition 3.1. For 71 ,..., 7 ^ G C with Sijj >0 (1 ^ j ^ r), the desingularized 

multiple zeta-function, which we also call the desingularization of (r((sj); (jj)), is 
defined by 


cr((^i);(7i)) 

(-ir 


n 


1 


• - lim , N I I / O • N / N I 

cgm\{i} 


^dtk (3.2) 


k=l 


for {sj) G C’’, where C is the Hankel contour used in (2.1). Note that (3.2) is well-defined 
because the convergence of the contour integral and of the limit with respect to c —1 can 
be justified from Theorem 3.4 (see below). 


Remark 3.2. By (2.1) and (1.14), we may say that equation (3.2) is a rigorous way to make 
sense of the nonsense equation (3.1). 


Example 3.3. In the case r = 1, set (r, 71 ) = (1,1) in (3.2). Similarly to |24l Theorem 4.2], 
we can easily see that 


C^(s;l) = lim 


(- 1 ) 


1 


1 c - 1 (e 2 "® - l)r(s) Jq Ve* - 1 e'’* - 1 


1 


^dt 
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m=l 


cm) 


= lim 


(- 1 ) 




(3.3) 

c->-l C — 1 ^ ’ 

Hence 1) can be analytically continued to C. As was mentioned in Introduction, this 

is the ’’proto-type” of desingularization. 

More generally we can prove the following theorem. 

Theorem 3.4. For 71,..., 7^ G C with iR'yj >0 (1 ^ j ^ r), 

C?“((sf);(7i)) 


n 


1 


U - l)r{sO Je-i™ 

1 


(7j); c) Yl tl^ ^dtk 
^ ^ k=i 


n 


(e 2 ^is, _ l)r(sk) 

(- 1 ) 


n lim 

C-5-1 




c 1 \ exp (jj tk) - 1 exp (cjj J2k=j 4 ) - 1 / k=i 




which can be analytically continued to C” as an entire function in {sj). 


For the proof of (3.4), it is enough to prove that if |c — 1| is sufficiently small, then there 


exists a function F : C” —>■ M>o independent of c such that 

|(c- l)-”^.((t,);(7,);c)| ^ F{{t,)) ((t,) E e”), 

f n K'°~^dtk\ < 00. 


(3.5) 

(3.6) 


k=l 


Now we aim to construct F{{tj)) which satisfies these conditions. Let (N’(e) = {z E C | |2;| ^ e} 
and §{0) = {z E C I | argz| ^ 6}. 

Let 71 ,... ,7r E C with IR'jj >0 (1 ^ j ^ r). Then the following lemma is obvious. 

Lemma 3.5. There exist e > 0 and Q < 9 < n/2 such that 

r 

7i^4 eN(1)US(0) (3.7) 

k=j 

for any (tj) E C”, where G is the Hankel contour involving a circle around the origin of radius 
e (see {2.1}). 

Further we prove the following lemma. 

Lemma 3.6. Let c E M \ {1} satisfying that |c — 1| is sufficiently small. Then there exists a 
constant A > 0 independent of c such that 

1 c 


c — 1 


-1 


ev - I e^y - I 


< 


Ae-^y/^ 


(3.8) 


for any y E N(l) U S(0). 
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Proof. It is noted that there exists a constant C > 0 such that 

1 c 


c — 1 


-1 


<C (yGK(l)), 


ev -I -I 

where we interpret this inequality for y = 0 as that for y —)• 0. Also, for any y G §(9) \ ^(1), 
we have 


|c — 1 


-1 


1 


eV - I e^y - I 


= |c — 1 


= |c — 1 


-1 


-1 


o^y _ npV 


cey + c — 1 


{ey - l){e'^y - 1) 
e^y - ey + {1 - c){ey - 1) 


{ey - l){e‘^y - 1 ) 


^ |c — 1 


-1 


cCy _ py\ 


1 


\ey-l\\e^y-l\ |e''2^-l| 


Hence it is necessary to estimate 


|c — 1 


-1 


^cy — py\ 


ley - iWe^y - 1 | 


We note that 


e^y - 1 


CXD j — j 


IE 


j! 


°° I 1 ^ 


Since |y| ^ JRy/ cos0, we have 


|c — 1 


-1 


<=ry _ pyl 


1=0 


1 


;(c-l)y _ ]^| 


|ey _ l||ecj/_ l| \i - e-y\\e^y - i\ \c - l\ 

1 






|1 - e-yWe'^y - 1 | 

|y|g%(|c-l|/cos6») 


|y|e 


\{r-i)y\ 


|1 - e-y\\e‘^y - 1| 

Therefore, if |c — 1| is sufficiently small, then there exists a constant A > 0 such that 


c — 1 


ley - llle^y - 1 


This completes the proof. 


Proof of Theorem 3.4. With the notation provided in Lemmas 3.5 and 3.6, we set 
F{{tj)) = A^fl exp I -SR(7 j ^ 4/2) | = exp [ - ^ ^ 4/2) 


i=i 


k=j 

k 


i=i 


k=j 


^'’exp -^3f?(4(^7j/2)) =A'’JJexp -3f?(4(^7^/2)) 


k=l 


i=i 


k=l 


i=i 


□ 


Then it is clear that F{{tj)) satisfies (3.5) and (3.6). Hence, by Lebesgue’s convergence 
theorem we see that (3.4) holds. 


Similarly to the proof of Theorem 2.1, since the contour integral on the right-hand side 
of (3.4) is holomorphic for all (s^) G C'’, we see that Cr^®((sj); (7^)) can be meromorphically 
continued to and its possible singularities are located on hyperplanes = 
r) outside of the region of convergence because — l)r(sfc) does not vanish at G Z^o- 
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Furthermore, for = Ik ^ N, the integrand of the contour integral with respect to tk on the 
right-hand side of (3.4) is holomorphic around tk = 0. Therefore, for Ik £ N, we see that 

1 c 


f (- 1 ) 

lim / lim- 

Jq c— >1 c 1 


exp 7i JZ=i tu) -1 exp cjj tu) - I 


lim 


'Ce 


e 11 gxp YJ- tj] - 1 exp (c 7 j Y]u=j “ 1 


t^k ^dtk 


t^k 


= 0 , 


because of the residue theorem, where = {ee*® | 0 ^ ^ 27r} for any sufficiently small 

e. Consequently this implies that ( 7 i)) has no singularity on Sk = h, namely 

( 7 j)) is entire. Thus we complete the proof of Theorem 


3.4 


Theorem 3.7. For 71 ,..., 7 ^ G C with iR'jj >0 (1 ^ j ^ r), 


n 

i=i 


1 - 7i Efc =7 exp 7 j tk) - I 


- 1 


(exp ( 7 ^ tk 


E (7,)) n ^ 


mi,...,mr=0 


-- rrii 


i=i 


Hence, for (kj) £ Ng, 

r 


Z=1 


E n ^i+Er=7 


^31 1 


j=\ 

"12. "22S*0 

I'lr, ..., izrr^O 




Proof. By (1.13), we have 
(- 1 ) 


lim , 

C- 7 I (c — 1 )’ 


■^r-((ti)Li;(7i)i=i;c) 


lim f[ Ell 

c-s-l J- J- C — 1 
i=i 


exp 7j EL, 4-1 exp cy^- ■ 4-1 


n 

i=i 


1 - 7i EL,' 4 exp 7 j Y]k=, 4-1 


exp 7j EL,' 4-1 


Hence we obtain (3.9) from (3.4). Also, by (1.13), we have 

(-L 


lim , 
c->-i (c — 1 ) 


:i3r((4)Li;(7i)Li;c) 


□ 


(3.9) 


(3.10) 
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( 

limTT y 

j=l ymj=l 


- 1 E[=j 

Bmj 


rrijl 


( 


n 

j=l \rnj=l 

( 




{ruj — 1 )! 


n 


OO ^1^33 Uj 




u^A UjA 




/ 


E ni'^ 

^11^0 i=i 
^12’ ^22^0 

^2r >•• •’ ^rr'^0 


i+EL, 7i 




Ed=i 




Hence, substituting the above relation into (3.4) and using the residue theorem with 

lim - 1 ) r(s) = {kGNo), 


s^—k 


we have 


c?“{(-%)47j))= n 

1=1 


X (2^*)" E n y 


i+EL, 7i 


EL, 1 


jl = l 

^ 12 ’ *^ 22^0 




Thus we obtain the assertion. 


□ 


Now we give an expression of ((sj); (,'lj)) in terms of Criisj); (1); ( 7 j)), which can be 
regarded as a multiple version of y®(s; 1) = (1 — s)C,{s) in the case r = 1 (see Examples 


3.3 


and 4.1). 


For Sj G C with 3?Sj > 1 (1 L i ^ r) and 71 ,..., 7 r G C with 3 ? 7 j > 0 (1 L J ^ r ), we set 

1 „ ” 


Ic,r ('Si, . . ■ , Sr j 7l) • • • ) 7r ) • — 


HdtjUt; 


nj=ir(sj) y[o,oo)’'^y yy 


n 


y V exp ( 7 j J2l=i 4-1 exp cjj J2l=i 4 ) - L 


(3.11) 


From Definition 3.1, we see that 




i (c - 1 )' 


r 4 ,r((sj); ( 7 i))- 
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For indeterminates Uj,Vj (1 ^ ^ r), we set 

r 

(Vj)) := n + • • • + UrVr)ivJ^ “ 


(3.12) 


i=i 


with the convention Uq ^ = 0 , and also define the set of integers {ai,Tn} by 


Si{Uj),{Vj)) = 


L- rrij 


(3.13) 






i=i 


where the sum on the right-hand side is obviously a finite sum. Note that the condition 
X)j=i ~ b Ibe summation indices m = ( rrij ) can be deduced from the fact that the 
right-hand side of (3.12) is a homogeneous polynomial of degree 0 in ( vj ), namely so is that 
of ( [3l3l ). 

Theorem 3.8. For 71 ,..., 7 ^ G C with iR'yj > 0 (1 ^ j ^ r), 

r 

Cr''''((Si); (7j)) = Y (^l,m{Y\{Sj)lj'jCr{si+mi,...,Sr + mr-,{l);{'yj)) (3.14) 


m = {mj ) 


i=i 


holds for all (sj) G C'", where (s)o = 1 and {s)k = s(s + 1) ■■■ (s + k — 1) (k G N) are the 
Pochhammer symbols. 


We emphasize here that each term of the right-hand side of (3.14) is meromorphic with 
infinitely many singularities but taking the above finite sum of the shifted functions causes 
‘miraculous’ cancellations of all the infinitely many singularities to conclude an entire func¬ 
tion. 


Remark 3.9. In (3.14), the condition ~ b implies that all zeta-functions appearing 

on the both sides have the same weight si 


Proof of Theorem \3.8[ First we assume that iftsj is sufficiently large for 1 ^ j ^ r. From 
( 2 . 1 ) with (f^j) = ( 1 ), we have 

.Sa —1 r 


Cr((sj);(l);(7i)) = 


1 


n 


*? 


0^=1 r(sj) 7(0, 00 )- exp(^7j Y.k=j “ 1 i=i 


Yldtj. 


(3.15) 


Using the relation 


lim-^^(^— 

c^i c — 1 — 1 


ecy - 1 
-1 + ey - yey 1 


ye^ 


we have 


{ey - 1)2 ey -1 (ey - 1)2 


(- 1 )^ 


= ^(y) (say), 


Cr^'^iiSj); ilj)) = lini - - :—Ic,r({Sj); (jj)) 

c^l [C — 1 ) 
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= lim ^ / TT dtj IT t 

(- 1 )/ 1 


Sj-l 

3 


n 

i=i 


c — 1 


exp 7i ELj 4-1 exp c-fj ^. 1=3 4-1 


(3.16) 




11^=1 r('Si) 7[o,oo)’- 


j=i k=j 


We calculate the last product of (|3.16[). Using the relations 


- 1 


= E= 

n=l 


-ny 


{ev - 1)2 


= 

n=l 


-ny 


we have, for J C {1,..., r}, 


n^jri'r n 


n 


7j 4 ) exp ( 7 , Yl.j 4 ) 


'|o,oor,.i ,=1 ' rfjexp(7,j;i,,4) - Dgj (exp(7jEU4)- 1 )' 

T V 00 T 

^[0,00) n(tJni = l k=j 

00 


HE rij exp^—nj 7 j ^ 4 ) ( 7 ^ tk 


j&J nj = l 


E (n 


njlj 


k=j j&J k=j 

r r r 


ni,...,nr71 JSJ 

3 i“l 

E (n(E nklk - ^ nfc7fe 

ni,...,nr71 JSJ k=l k=l 


n dtj ]J f/ ^ Yl exp (-4 ^ nk-fk) n {Y 1 

4~)’'j=i j=i j=i ■ - ■ 


k=l j&J k=j 


r r r j r 

n dtj exp (-4 Hklk) n {Y1 


Jio,ocrj^^ j=i j=i 

3 i~i 

E'-^ E (n(E nklk - Y 

iSNg ni,...,nr71 iSJ k=l k=l 


J n^4n4''^'' 'fiexp(-4^ 

i[0,oo) 


k=l j&J k=j 


nklk 


j=i j=i j=i 

3 i-i 

'nklk 

^SNq ni,...,n^71 j£J k=l k=l 


3 3 -^ r 

(n(z] nklk - Y ))nr(e 7 + i 7 ) 


1 


i=i 


(s 


E rirfe + '^l E E - 

im^ j=l ni,...,nrllKcJ\{l} j=l {Y.k=l'^klk 


k=l nklk 

1 


^j+l~^jeJ\K-^j+ieK 


I 

Yhi Y + + (3-17) 

ieNg KCJ\{1} j=i 
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where \K\ implies the number of elements of K, 

'l {iGl) 


— 


0 


for I C J, and 




(3.18) 


j&j k=j j=i 


Hence, by (3.16) we have 


E E E (- 1 )'^' (ri “ <^i+ieA:); ( 1 ); (7i)) 




l£% KCJ\{1} 


./C{l,...,r} ACJ\{1} leNg j=l 

Finally we set 


(3.18) 


-.= Y. E EUill 

JC{l,...,r} ACJ\{1} I&Wq i=l 

and aim to prove that 

S((«i), = H{{uj), {vj)). 


„,h„,h-SjeJ\K-Sj+ieK 


(3.20) 


It follows from (3.18) that 


H{{uj),ivj))= Y 1 E (-i)'^^'^'(nE 


-^jeJ\K-^j+ieK 


Jc{l,...,r} KCJ\{1} 

JC{l,...,r} j&J k=j 


jeJ k=j j=l 


E n (-^ 1 ') n ^A- 

KcJ\{l}j&J\K j£K 


Since Uq ^ = 0, we have 


E n (-E)n"7-i=n(-7‘+v7 

^cJ\{i}ieJ\A: jeK j&j 


Hence we obtain 


(vj)) = E n (E ^ 

JC{l,...,r}ieJ k=j 


n ((^ (-^^/+v'l)+ 1 ) 

j=l k=j 
r r 

n(^ “ {^'^kVk){vJ^ - v\)) = 9i{Uj), (Vj)). 


(3.21) 


i=i 


k=j 


Combining (3.13), (3.19) and ( 3.21[ ), and regarding {sj)i- and Cr{{sj + Ij)', (1); (7j)) as inde- 
terminates and v[^, we see that (|3.14|) holds when 3f?s,' is sufficiently large for 1 ^ ^ r. It 

is known that each function on the right-hand side can be continued meromorphically to C'’ 
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(see [m Theorem 1]). Since i'lj)) is entire, we see that (3.14) holds for all (sj) G C^. 

□ 


Thus we complete the proof of Theorem |3.8[ 

4 . Examples 

Example 4.1. In the case r = 1 and 71 = 1, we have 

9{ui,Vl) = 1 - UlVlV^^ = 1 - ui, 
namely, ao,o(l) = 1 and ai,o(l) = —1- Hence we have 


1 ) = lim = (s)oCi('S; 1; 1) - ('S)iCi('S; 1; 1) = (1 - s)C(s), 


(- 1 ) 


c^l C — 1 


which coincides with (3.3). We see that 

Cf-(1;1) = -1 

and 

C^(-fe;l) = (-l)'=Hfc+i (A:GNo). 

Example 4.2. In the case r = 2, we can easily check that 

g((Uj), (uj)) = (1 - {uiVl + U2V2)v^^){l - U2V2 {v2^ - V^^)) 
= (1 — ui)(l — U2) + (^2 — UlU2)v^^V2 — 


Then (3.14) implies that 

C2®"('Si,'S2;7i, 72) = (Si - 1)(S2 “ l)C2(si,S2; (1);71,72) 

+ S 2 {S 2 + 1 - si)C2(si - 1 ,S 2 + 1 ; ( 1 ); 7 i, 72 ) 
- S 2 {S 2 + l)C 2 (si - 2, S2 + 2; ( 1 ); 7 i, 72 ). 


Let k,l G Nq. By (3.10), we obtain 


C2"'^«(-A:,-/;7i,72) = (- 1 )"+^ ^ 

Remark 4.3. Setting ( 71 , 72 ) = (1,1) in ( |4.1[ ), we obtain 
C2'''(si, S 2 ; 1,1) = (si - l)(s2 - l)C2(si, S 2 ) 

+ 82(^2 + 1 — 'Si)C 2 ('Sl — 1 , S 2 + 1 ) ~ S2{S2 + l)C 2 ('Sl — 2 , S 2 + 2 ). 

From Theorem 


(4.1) 


(4.2) 


(4.3) 


3.4 


we see that C 2 (si, S 2 ; 1,1) on the left-hand side of ( |4.3| ) is entire, though 
each double zeta-function (dehned by on the right-hand side of ( |4.3[ ) has infinitely 

many singularities (see (0.4)). In fact, we can explicitly write —n; 1,1) in terms of 

Bernoulli numbers by (4.2), though the values of C 2 (si,S 2 ) at non-positive integers (except 
for regular points) cannot be determined uniquely because they are irregular singularities 
(see d]). 

Example 4.4. In the case r = 3, we can see that 

g((uj), (vj)) = (1 - {uivi -h U2V2 + U3i;3)uf ^)(1 - {U2V2 + U3'i;3)(u^^ - uf ^)) 


X (1 - U3V3{v^ ^ - V 2 ^)) 


.-l^ 
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= -{Ui - 1)(U2 - l)(lt3 - 1) + (itl - 1){U2 - m)u3V2^V3 
+ (mi - l)ulv2^vl + (iti - U2)U2{U3 - l)v'{^V2 
+ U3{-Ui + 2u2 - U1U2 + U2 + U1U3 - 2u2U3)v^^V3 
— u\{ — l + Ui — 2 u2 + U3)v^^V2^V3 + 

+ ul{u3 - l)v^‘^vl - U2{2 + U2- 2 u3)u3V^‘^V2V3 
+ u\{ — l — 2 u2 + U3)v^‘^V 3 — ulvi‘^V2^V3. 


Therefore we obtain 


C3''''('S1>'S2,S3;7U 72,73) 

= -(si - 1)(S2 - 1)(S3 - l)C3(si,S2,S3;(l);7l,72,73) 

+ (si - 1)(-1 + S2 - S3)S3C3(S1,'S2 - 1, S 3 + 1; (1); 7l> 72,73) 

+ (si - 1 ) 53(53 + 1)C3(5i,52 - 2 ,S3 + 2; (1); 71 , 72 , 73 ) 

+ (-1 + Si - S 2 )S 2 (S 3 - 1)C3(5i - 1,S2 + 1, S 3 ; (1); 7l, 72,73) 

+ S 3 (S 2 - S 1 S 2 + sl + S 1 S 3 - 2 s 2 S 3 )C 3 ( 5 i - 1, S 2 , S 3 + 1; (1); 7l, 72,73) 

- 53(s 3 + 1)(1 + Si - 2 s2 + S3)C3(si - 1,S2 - 1, S3 + 2; (1); 7l, 72,73) 
+ 53 (s3 + 1)(S3 + 2)C3 (si - 1 ,S2 - 2, S3 + 3; (1); 71,72,73) 

+ S 2 (S 2 + 1)(S3 - l)C3(5i - 2,S2 + 2, S 3 ; (1); 7l>72,73) 

- S 2 (l + S 2 - 2 s 3 )s 3 C 3 (si - 2, S 2 + 1, S 3 + 1; (1); 71 , 72 , 73 ) 

+ 53(s3 + 1)(1 - 2s2 + 53)C3(5i - 2, S 2 , S 3 + 2; (1); 7l, 72,73) 

- S 3 (S 3 + 1 )(S 3 + 2 )C 3 (si - 2,S2 - 1, S 3 + 3; (1); 71 , 72 , 73 ). 


Let A:, Z, m G Nq. By (3.10), we have 

^des 


m m—u I 

CS“(-fc,-!,-m;7l.72,73) = E E 

i/=0 p=0 K,=0 


m 


Kj \v p 


sy TD TD D I —t^-\-p+'^ Tn — V — p+\ 


where (™) = , ,, 

\v p) u\ p\{ra—v—p)\ 


-p)\ 


Remark 4.5. Our desingularization method in this paper is for multiple zeta-functions of the 
generalized Euler-Zagier type (0.1). In [3, we will extend our desingularization method to 
more general multiple series. 


Remark 4.6. Arakawa and Kaneko [3] defined an entire function ^fc(s) associated with poly- 
Bernoulli numbers mentioned in Remark 


1.12 


It is known that, for example. 


6 ( 5 ) = sC(s + 1), 

? 2 (s) = —C 2 ( 2 , s) + sC 2 (l, s -|- 1 ) -|- C(2)C(s). 

Comparing these formulas with (|3.3|) and (|4.1|), and using the well-known formula 


C 2 ( 0 ,s)= 


1 




N=2 
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we obtain 

and 

(1 - 5 ) 6 ( 5 ) = C 2 ^'=^( 2 , 5 ; 1,1) - cf ^(2)C^(5; 1) - (5 + l)Ci"^^(5 + !;!) + sC^(5 + 2; 1), 
C2^'^^(2, 5; 1 , 1 ) = (1 - 5)6(5) + 6 (i) 6(5 - 1 ) - (5 +1)6(5) + 56(5 + 1 ). 

Note that the both sides of the above relations are entire. It seems quite interesting if we 
acquire explicit relations between 6 ( 5 ) and Cfc'^®(( 5 j); (1)) for any k ^ 3. 

Related to the Connes-Kreimer renormalization procedure in quantum field theory, Guo 
and Zhang m and Manchon and Paycha |16) introduced methods using certain Hopf algebras 
to give well-defined special values of the multiple zeta-functions at non-positive integers. 


Example 4.7. According to their computation table (in loc.cit.), Guo-Zhang’s renormalized 
value C^^(O) “2) of 6 ( 51 , 52 ) at its singularity (si, 52 ) = (0, —2) is 


1 


C 2 (0,-2) =-, 

S 2 V > ; ^20’ 


while Manchon-Paycha’s value C 2 "(0, —2) is 


cr(o,-2) = 


720' 


On the other hand, our desingularized method gives 


C 2 "-( 0 ,- 2 ; 1 , 1 ) = -, 

so these three methods give values different from each other. 


Question 4.8. Are there any relationships between our desingularization method and their 
renormalization methods? 


Finally we emphasize that since our ^)?®®((sj); (1)) is entire, their special values at integer 
points which are neither all positive nor all non-positive are well-determined. These values 
might be also worthy to study. We conclude this paper with the announcement of explicit 
examples of those values. 


Example 4.9. We have 

CfVl.4;l.l) = C(3)-C(4), 

Cf(3,-3il.l) = 65|«=*)' 

C24=-(4,-3;1,1) = 6^«2)-Pc(4). 

Also we can give the following examples for non-admissible indices: 

C24“(2,1;1,1) = -C(2) + 2C{3), 
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Cf^(3,l;l,l) = 2C(3)-^C(4). 

For the details, see [9]. 
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